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There exists a growing evidence for the anomalous transparency of the Universe for energetic 
gamma rays. Popular explanations include conversion of photons into hypothetical axion-like parti¬ 
cles (ALPs) and back in astrophysical magnetic fields. Two distinctive scenarios of this conversion 
have been put forward: either it happens in the (host galaxy of the) gamma-ray source and in the 
Milky Way, or the photon-ALP oscillations take place in the intergalactic magnetic fields all along 
the way between the source and the observer. Here we point out that, given recent astrophysical 
constraints on ALPs and on intergalactic magnetic fields, these two mechanisms imply very different 
ALP masses and couplings. Therefore, confirmation of the anomalies and identification of one of 
the scenarios would mean cornering of ALP parameters to a particular narrow region. We discuss 
approaches to distinguish between the two mechanisms and present some preliminary indications in 
favour of the galactic scenario. 
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I. INTRODUCTION 

Since the early years of very-high-energy (VHE, above 
100 GeV) gamma-ray astronomy, more and more distant 
sources have been discovered in this energy band, see e.g. 
Refs. [IH3]. This was surprising because gamma rays of 
these energies should interact with the background ra¬ 
diation to produce electron-positron pairs [1], the pro¬ 
cess which results in a strong suppression of the flux 
from distant sources. The apparent weakness of the 
suppression might be attributed either to overestima¬ 
tion of the amount of background radiation or to pecu¬ 
liar emission or absorption mechanisms at work in par¬ 
ticular sources. However, recent studies indicate that 
these explanations are hardly relevant: modern models 
of extragalactic background light, e.g. Refs. 0 0: used 
in the studies, already saturate lower bounds [7] from 
simple galaxy counts; while analyses of sets of distant 
sources [HHin] revealed unphysical redshift dependence 
in the required peculiar features of the emission spectra. 

One of the most intriguing explanations of the ensem¬ 
ble of the data is that we might observe effects of a new 
hypothetical axion-like particle (ALP) (see Ref. [TT] for 
a review). Depending on parameters of the ALP and on 
values of the intergalactic magnetic field (IGMF), two 
scenarios may work. Photon-ALP oscillations in IGMF 
may effectively increase the photon mean free path Hi- 
Ill]. Alternatively, if IGMF is weak, photons may effi¬ 
ciently convert to ALPs in the source [TS] or in its neigh¬ 
bourhood, cluster or filament [16] , with reconversion back 
to photons in, or close to, the Milky Way. The main pur¬ 
pose of this study is to emphasise differences between the 
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two scenarios and to demonstrate that, given recent as¬ 
trophysical constraints, they correspond to distinguish¬ 
able ALP parameters. The key new ingredients con¬ 
tributing here are the reanalysis of the supernova (SN) 
1987A data [T7| and new upper limits on intergalactic 
magnetic fields [TH|. We also present some evidence in 
favour of the galactic scenario and discuss future obser¬ 
vations which may help to discriminate between the two 
options (see also Refs. mm)- 
The rest of the paper is organized as follows. In Sec.jlll 
we briefly review the evidence for the anomalous trans¬ 
parency of the Universe for gamma rays. In Sec. m 
we describe two scenarios of the ALP-photon mixing 
which may explain the observed evidence. Then, we 
compare required ALP parameters, calculated by means 
of detailed analyzes elsewhere, with recent astrophysical 
bounds and present argumentation for the main point of 
the paper. In Sec. |IV[ we consider ways to distinguish 
the two scenarios and present some evidence in favour 
of one of them. Future prospects to test the ALP ex¬ 
planation of the anomalous transparency of the Universe 
for gamma rays and to single out one of the scenarios, 
thus constraining the ALP parameters, are discussed in 
Sec. El 


II. EVIDENCE FOR ANOMALOUS 
TRANSPARENCY 

The modern evidence for the anomalous transparency 
of the Universe for energetic gamma rays is based on 
studies of ensembles of distant VHE sources. The ob¬ 
served spectra of these sources have been corrected for 
the pair-production effects (“deabsorbed”) within the 
most conservative, i.e. lowest-absorption, models, to ob¬ 
tain intrinsic spectra emitted at the sources. These in- 
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trinsic spectra exhibit unphysical redshift dependence 
which is readily interpreted as an overestimation of the 
absorption even in the minimal models. 

Upward breaks, or unusual spectral hardenings, have 
been found in deabsorbed spectra of many individual 
sources (see e.g. Ref. 0)- Statistically significant hard¬ 
ening never presents in observed spectra which, contrary, 
often exhibit mild softening at high energies. In 2012, 
Horns and Meyer analysed [8] a sample of 7 blazars ob¬ 
served at optical depths r > 2 with respect to the pair 
production. The blazars were observed by imaging at¬ 
mospheric Cerenkov telescopes (lACTs) and have red- 
shifts z < 0.536. They found an evidence that positions 
of the upward breaks in deabsorbed spectra of blazars 
change with redshift in such a way that they always occur 
at the energy where the absorption becomes important. 
This was surprising because astrophysical properties of 
blazars in the sample were very similar for close and dis¬ 
tant sources. The probability that this effect is a chance 
fluctuation, estimated in Ref. [8] by a statistical proce¬ 
dure based on the Kolmogorov-Smirnov test, corresponds 
to that of a 4.2cr fluctuation in a Gaussian distribution. 

Rubtsov and Troitsky [9] considered in 2014 a sample 
of 20 blazars observed at optical depths r > 1. lACT 
results have been supplemented with the FERMI-LAT 
data, which allowed to extend the redshift range up to 
z « 2.156. Assuming the breaks in deabsorbed spectra 
at energies for which r = 1, we found that the strength of 
the break, that is the difference between the power-law 
spectral indices below and above the break point, is a 
function of the redshift and does not depend on the phys¬ 
ical properties of a blazar. The statistical significance of 
this unphysical dependence, indicating overestimation of 
the absorption and therefore an anomaly, has been cal¬ 
culated in Ref. [9] by means of the usual analysis and 
corresponds to a 12Aa fluctuation of a Gaussian. 

In 2015, Galanti et al. [10] considered a sample of 39 
blazars {z < 0.536) detected in VHE gamma rays, inde¬ 
pendently of the opacities tested. They described deab¬ 
sorbed spectra as power-law functions and did not fit the 
spectral breaks. The power-law spectral index was found 
to be redshift-dependent, which is not expected in any as¬ 
trophysical model and again indicates the anomaly. This 
result confirms the observation of Ref. jH] though the in¬ 
trinsic scatter of spectral indices and the limited redshift 
range make the result less pronounced than in the break 
study. 

The significance estimates quoted above are based on 
statistical analyses only and are subject to systematic 
uncertainties, which are discussed in Refs. 00. In par¬ 
ticular, Ref. 0 shows that under the worst assumptions 
about systematic errors, the significance of the effect is 
reduced by ^ 1.6 ct. Ref. 0 presents several tests of the 
robustness of the result; however, a detailed quantitative 
study of systematic uncertainties is hardly possible with¬ 
out a complete sample of sources (tests of the Malmquist 
bias have been presented in Ref. 0 for the part of the 
sample taken from the FERMI-LAT catalog). 


Recently, Biteau and Williams jUj criticised the re¬ 
sults of Horns and Meyer 0 and claimed they do not 
see any anomaly in the ensemble of deabsorbed spec¬ 
tra of 38 blazars observed by lACTs. Since they did 
not select the blazars by the energies at which the data 
are available, their sample is dominated by the sources 
for which the absorption is, and should be, low or neg¬ 
ligible (as we know from Ref. 0 , only 15 blazars were 
firmly observed by lACTs at r > 1, and 0 only 7 at 
T > 2). In addition, they used multiple spectra for a 
given blazar in their work; this introduces additional 
statistical weight to better-studied nearby, and there¬ 
fore unabsorbed, sources. Even for heavily absorbed 
sources, they apply a different method to derive their own 
model of background radiation and the intrinsic spec¬ 
trum. All these points might explain the discrepancy 
between Refs. [21] and 0 which Biteau and Williams 
attribute to the use of the Kolmogorov-Smirnov test by 
Horns and Meyer. In any case, this critique is irrele¬ 
vant to our work 0 in which we made use of the usual 
test to determine the significance of the anomalous 
transparency in a clean sample. 

The only astrophysical explanation of these anoma¬ 
lies |22] requires some non-conventional assumptions. It 
assumes that a sufficient number of ultra-high-energy cos¬ 
mic protons are accelerated in precisely the same gamma- 
ray blazars. Unless extragalactic magnetic fields are as 
low as < 10“^^ G everywhere along the line of sight, 
this scenario may have tensions with the observation of 
fast variability of 4C-I-21.35 at very high energies [23] . 
In what follows, we assume the anomalies are real and 
concentrate on a different kind of their explanations. 


III. EXPLANATIONS 

The problem of unphysical spectral breaks may be alle¬ 
viated if the gamma-ray attenuation is reduced by means 
of some mechanism. Indeed, the breaks are seen precisely 
at the energies where the correction for attenuation be¬ 
comes important (lower energies for larger distances); re¬ 
duction of the attenuation diminished the correction and 
removes the breaks 0 : given observed softening of high- 
energy spectra, one may even obtain expected power-law 
shapes if the absorption is present but reduced. However, 
the usual deabsorption procedure is based on firm stan¬ 
dard physics and very conservative assumptions about 
the photon background, so only new-physics effects might 
help. Besides the possibility of the Lorentz-invariance vi¬ 
olation, the only known explanation involves ALPs. 

An ALP mixes with photons in external magnetic 
field |24j . which may allow to suppress the attenuation 
due to pair production: gamma-ray photons convert to 
ALPs, then travel unattenuated and eventually convert 
back to photons. The photon beam is still attenuated, 
but the flux suppression becomes less severe. A useful 
collection of formulae describing the mixing for astro- 
physical environments, as well as a quantitative discus- 


3 


sion of various scenarios, may be found in Ref. [16) . To 
reduce the opacity of the Universe for TeV gamma rays 
from blazars, two particular scenarios involving ALPs are 
important. The purpose of the present study is to em¬ 
phasise and to explore the difference between the two 
approaches. 

The first scenario implies that the intergalactic mag¬ 
netic field is strong enough to provide for ALP/photon 
conversion all along the path between the source and 
the observer. Originally suggested in Ref. m in a 
different context, this mechanism, known also as the 
DARMA scenario, was invoked for the TeV blazar spec¬ 
tra in Ref. m- If it is at work, then the photon/ALP 
mixed beam propagates through the Universe and, since 
the photons are attenuated while ALPs are not, the ef¬ 
fective suppression of the flux is smaller compared to 
the pure-photon case. It is easy to demonstrate that, 
for a sufficiently long propagation through domains of 
randomly oriented magnetic fields, the optical depth is 
effectively reduced by 2/3 in this scenario. A more de¬ 
tailed study is given, for instance, in Ref uni, from where 
the most recent constraints on the relevant ALP pa¬ 
rameters are obtained: the ALP mass m < 10“® eV 
and the ALP-photon coupling ga-yj determined from 
^ = {B/nG){ga-y-y x 10^^ GeV) > 0.3, that is ga-y-y // 
3 X 10“^^ GeV“^. In what follows, we will refer to this 
mechanism as the “intergalactic conversion” and use the 
parameter constraints |I0j for this scenario. 

The second approach assumes that there are quite 
strong magnetic fields inside or around the source, as 
well as around the observer, while for the most part of the 
distance the beam travels in weak magnetic fields, insuffi¬ 
cient for ALP/photon mixing. The conversion may hap¬ 
pen either in the blazar itself and in the Milky Way [15] 
or in the galaxy cluster or filament m ( see also a more 
detailed subsequent study in Ref. [IS]) containing the 
source and the observer, in various combinations. To get 
a qualitative idea of the effect of this mechanism on the 
gamma-ray attenuation, one might consider the case of 
maximal mixing, which is certainly an unrealistic over¬ 
simplification for many particular sources. Then 1/3 of 
the original photon flux is converted into ALPs close to 
the source while the remaining 2/3 attenuate in a usual 
way. Close to the observer, 2/3 of the ALPs convert back 
to gamma rays and may be detected. A more detailed 
recent study of this mechanism is presented in Ref. m, 
where it is called “the general-source” scenario. Notably, 
this scenario requires ga-y-y > 2x 10“^^ GeV“^ because for 
lower values of the coupling, the path of the ALP-photon 
beam would be too short for efficient conversion even for 
maximal mixing (the Hillas-like argument). In the rest 
of the paper, we refer to this mechanism as the “galac¬ 
tic conversion” and use parameter constraints derived in 
Ref. [in] for this cas^ 


^ In Fig. [^ below, the corresponding line was smoothed with re¬ 
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FIG. 1. The ALP parameter space (ALP-photon coupling 
ga-y-f versus ALP mass m) with current constraints (see text). 
Regions corresponding to the Galactic m and intergalac¬ 
tic [12] ALP/photon conversion explanations of the gamma- 
ray anomalies are indicated; they extend to the excluded re¬ 
gions as shown by arrows. Given all constraints, the two 
regions are well separated. 

Clearly, if we knew the values of IGMF, it would be 
quite easy to choose between the two scenarios. Unfortu¬ 
nately, the values are a subject of strong debates, and we 
have to seek other ways to disentangle the two options. 

Parameters of the ALP, that is m and gayy, required 
for efficient operation of one or another mechanism, have 
been discussed in Refs. [MIIIIIS]; we refer the reader 
to these works for details of the derivation. Theoreti¬ 
cally, they overlap in a large range of allowed parame¬ 
ters. However, when the most recent experimental and 
astrophysical constraints are taken into account, the pa¬ 
rameter regions allowed for the two scenarios become dis¬ 
connected; this means that if we determine that one or 
another scenario works in Nature, we strongly constrain 
the ALP mass and coupling! We illustrate this fact in 
Fig.[T] where shaded blue areas, excluded by constraints 
from the CERN axion solar telescope (CAST, Ref. [26)1. 
evolution of the horizonthal-branch (HB) stars [57], re¬ 
analysis of the SN 1987A data m and HESS constraints 
from the absence of irregularities in a blazar spectrum 
[20) . indicate the most restrictive relevant limits. 

One may see that the key constraint contributing to 
the separation of the two regions is that of Ref. [Hi¬ 
lt was obtained from non-observation of the gamma-ray 
flare corresponding to the SN 1987A explosion by the 
Gamma-Ray Spectrometer of the Solar Maximum Mis¬ 
sion. The fluence limits used for the constraint are 99% 
CL; however, the responce of the instrument to pho¬ 
tons at large angles to viewing directions, relevant for 


spect to that of Fig. 4 of Ref. m since we believe the features 
presented there are overprecise, given observational uncertainty 
in GMF. 
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SN 1987A, is quite uncertain. The authors of Ref. m 
use results of the previously published Monte-Carlo sim¬ 
ulations. One may ask then, maybe the limit of Ref. m 
is wrong and the point of the present paper is destroyed? 
The answer is no. Indeed, for m < 10“® eV, where 
the intergalactic scenario may work with B satisfying 
the observational limits, the constraint m reads as 
Qa-y-y ^ 5 X 10“^^ GeV“^, while the galactic scenario re¬ 
quires gay-y ^ 2 X 10“^^ GeV“^, that is four times higher. 
Since the gayy limit changes as the fourth root of the flu- 
ence, the two inequalities for gayy might be brought into 
agreement by a ^ 4^ = 256 times error in the effective 
area. Clearly, even a rough estimate of the effective area 
of a satellite experiment can hardly be a factor of 256 
wrong. 

The separation of the two regions, which are often uni¬ 
fied in a single large band referred to as the “transparency 
hint” in relevant plots, is remarkable. It is instructive to 
compare this result with those of Ref. m which, if con¬ 
sidered superficially, suggest that the two mechanisms 
may work at the same time. However (see Fig. 4 of 
Ref. [in]), this might happen only in the “optimistic” 
scenario with the IGMF strength of 5 nG, now excluded 
by the study of Faraday rotations of distant sources m- 
We use parameters from the most recent study of the 
intergalactic scenario nni which, though published be¬ 
fore Ref. [im, does not assume IGMF in excess of the 
new limit, 1.2 nG. Hence, the probability of conversion 
in IGMF is lower and the parameter region is shrinked, 
compared to Ref. m- 

There exist some concerns about the possibility of ef¬ 
ficient axion-photon conversion in blazars, see Ref. |28j . 
The key ingredient in this reasoning is the Quantum Elec¬ 
trodynamics (QED) strong-held term which becomes im¬ 
portant at 

(i) (looliv) ~ <‘> 

see e.g. Ref. [TB] (the ALP-photon mixing in this regime 
may also be affected by photon-photon dispersion, see 
Ref. [29]). The concerns are however not critical for 
the galactic conversion scenario for the following rea¬ 
sons. Ref. [55] considers separately BL Lac type objects 
(BLLs) and flat-spectrum radio quasars (FSRQs). For 
BLLs, assuming magnetic helds B ^ (0.1 — 1) G at the 
site where gamma rays are produced, the condition 0 
is approached at the energies of interest, so the photon- 
ALP conversion probability may depend strongly on un¬ 
known details of the environment and becomes hardly 
predictable. However, if ALP parameters allow for the 
conversion in the Milky Way, one expects that the conver¬ 
sion in a BLL host galaxy or cluster is also allowecQ and 
the galactic scenario may work. The situation is more 


^ One should note, however, that the host galaxies of blazars are 
elliptical while the Milky Way is spiral. While the strength of 


contrived for FSRQs, for which, in Ref. [55], magnetic 
field values i? ~ (1 — 10) G are assumed (for the gamma- 
ray emitting region). The QED term starts to suppress 
the photon/ALP conversion at energies if > 10 GeV in 
this region, but efficient conversion is possible in outer 
parts of the galaxy and the cluster. It may be tricky 
for gamma rays above if ~ 20 GeV to get there because 
of intense pair production in the broad-line region, pre¬ 
sumably located between the gamma-ray emitting site 
and the lobe. However, one may point out that gamma 
rays up to several hundred GeV have been observed from 
ESRQs [TJ[53] (see Ref. [35] for a discussion of ALP expla¬ 
nation of the observation [53]), which means they escape 
the broad-line region somehow. Once escaped, thay may 
equally well convert to ALPs outside. 

IV. DISCRIMINATION BETWEEN GALACTIC 
AND INTERGALACTIC SCENARIOS 

Anisotropy. The magnetic field of the Milky Way 
galaxy has a complicated structure, and the probability 
of the ALP/photon conversion there, which is required 
in the galactic scenario, depends strongly on the direc¬ 
tion. Evidence for direction dependence in the anomalous 
transparency of the Universe may therefore be a strong 
argument in favour of the galactic scenario [T3[[I31I33]. 

In Ref. [15], it was pointed out that the positions of 
a few TeV blazars with redshifts 2 > 0.1 known by that 
time fit surprisingly well the regions in the sky where 
the conversion probability, calculated within the model of 
the Galactic magnetic field (GMF) of Ref. [34], are high. 
Here, we assume this as a hypothesis and attempt to test 
it with new observational data. Clearly, more elaborated 
approaches should be used in further studies (see Ref. [33] 
for an attempt on which we comment below and Ref. [35] 
for a different approach). We consider a sample of blazars 
with firm detection beyond r = 1 which consists of 15 
objects observed by lACTs and 5 objects observed by 
EERMI LAT (the sample of Ref. [9]), supplemented by 
additional 6 blazars rejected in Ref. [S] because of the 
insufficient number (4 with 5 required) of data points for 
fitting spectra with breaks. We drop 4 nearby objects 
with z < 0.1 from the sample, like it was done in Ref. [15j . 

Figure [^ represents the distribution of these objects 
in the sky together with the conversion probability for 
the same magnetic-field model [34]. Though there exist 
more elaborated modern GMF models [331137], the test 
of the original claim should be performed with the same 


the magnetic field in the Milky Way is assumed to be typical, 
elliptical galaxies contain fewer cosmic-ray electrons, and hence 
less synchrotron emission, making detection of the galaxy-scale 
magnetic fields tricky and, for today, uncertain m- Qualitative 
theoretical arguments suggest m small correlation length for 
the turbulent field. In any case, giant elliptical galaxies often 
reside in groups or clusters whose magnetic field also allows for 
the photon-ALP conversion. 
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FIG. 2. The skymap (Galactic coordinates, Hammer projec¬ 
tion) with positions of blazars with detected gamma-ray flux 
at energies for which r > 1 (red, 0.1 < z < 1; blue, 2 ; > 1), see 
text. Deeper shading corresponds to higher ALP-photon con¬ 
version probability in the Milky Way (the Galactic magnetic 
field model of Ref. |34p. 



^conv 


FIG. 3. The distribution of the values of the ALP-photon 
conversion probability in the Milky Way calculated for direc¬ 
tions to blazars with detected gamma-ray flux at energies for 
which T > 1 (red line) and for random isotropic directions 
(gray shading). 


one. The objects indeed follow the regions of high conver¬ 
sion probability, qualitatively confirming the trend seen 
in Ref. [TS]. Two of 22 objects are seen in the regions of 
low conversion probability; this may be well understood 
given uncertainties in the GMF models. 

It is not possible, however, to rigorously test the hy¬ 
pothesis quantitatively, because the blazars we discuss do 
not form a complete isotropic sample. While 5 FERMI- 
LAT sources were selected from a more or less uniform 
full-sky data, the remaining 17 objects were arbitrarily 
chosen for observations with pointed runs of small-field- 
of-view lACTs. Nevertheless, for illustration, we present 
here the results of a simple statistical test of the hypoth¬ 
esis, keeping in mind its qualitative level. For each of the 
22 sources in the sample, we calculate the ALP/photon 
conversion probability in the GMF of Ref. [34]. The same 
distributions were calculated and averaged for 100 sets of 
22 objects distributed isotropically in the sky. Figure]^ 
compares the distributions. The Kolmogorov-Smirnov 


probability that the distribution seen for the real data 
is a fluctuation of that for simulated directions is 0.02. 
We do not assign any statistical significance to this re¬ 
sult for the reasons mentioned above (incomplete sam¬ 
ple), though the entire picture does not contradict the 
galactic conversion scenario. 

Another approach to test the anisotropy effect was 
suggested by Wouters and Brun |33|. They proposed 
to study autocorrelation patterns in the directional dis¬ 
tribution of differences between FERMI-LAT and TeV 
spectral indices of blazars. We note that these differences 
are not equal to the spectral breaks discussed in Sec. |n| 
since the breaks happen at distance-dependent energies, 
not always between the bands considered in Ref. [33] . 
The procedure is motivated by the fact that, due to our 
limited knowledge of the GME, actual conversion prob¬ 
abilities vary strongly from one GME model to another, 
while all models predict patchy large-scale anisotropy in 
the distribution of these probabilities. However, as they 
show themselves |33| , the patterns in the autocorrelation 
function also vary strongly from model to model. Glearly, 
quantitative estimates of statistical significance would re¬ 
quire a complete sample of sources in this approach, sim¬ 
ilarly to a more direct approach m we discuss here. 

Distant objects. In the ideal case and in the long¬ 
distance limit, the effective optical depth talp behaves 
differently in the two scenarios: for intergalactic conver¬ 
sion, Talp (2/3)r (and therefore grows approximately 
linearly with distance, like the standard optical depth r), 
while for the galactic scenario, assuming maximal mixing, 
it reaches a constant, distance-independent value corre¬ 
sponding to the flux suppression by a factor of ^ 2/9, 
that is Talp ~ 1-5. At a certain redshift Zcrit, which 
value depends on the details of the absorption model 
and of magnetic fields assumed, the two suppression fac¬ 
tors are equal, while beyond 2crit, the effective absorp¬ 
tion becomes stronger and stronger in the intergalactic 
scenario, remaining essentially constant in the galactic 
one. This means that for very high redshifts, the anoma¬ 
lous transparency effect would hardly be seen in obser¬ 
vations for the intergalactic scenario, therefore any ev¬ 
idence for the effect for very distant sources speaks in 
favour of the galactic conversion. The blazars observed 
by lACTs have measured spectroscopic redshifts up to 
2 = 0.536 (3G 279; a lower limit of z > 0.6 exists for 
PKS 1424-1-240), and all anomalous-transparency effects 
derived from them are equally well described by both 
mechanisms. However, inclusion of much more distant 
blazars, for which the absorption becomes significant at 
energies in the FERMI-LAT band, changes the picture: 
the 12-sigma anomaly reduces to ^ 5cr (and therefore re¬ 
mains present) when intergalactic conversion is assumed, 
in addition to the usual absorption, but diminishes to 
~ 2(7 (and therefore disappears) for the assumption of 
the galactic scenario |9]. While these results have been 
obtained in a simplified approach, the difference between 
scenarios is so pronounced that it could hardly be re¬ 
moved by any detalisation of the analyzis. Prospects for 
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observations of distant gamma-ray sources are discussed, 
in the ALP context, in Ref. [38] , 

Intergalactic magnetic fields. The intergalactic 
scenario requires rather high IGMF, B ^ — 

10“®) G, otherwise the conversion probability would be 
too low. The value of IGMF is irrelevant for the galac¬ 
tic scenario provided the SN 1987A constraints on ALP 
parameters are satisfied. Present-day knowledge does 
not allow for a firm conclusion about the real values 
of R. A number of constraints are summarized in the 
review |39| . The most stringent observational limit, 
based on the redshift independence of the Faraday ro¬ 
tation from distant sources, is B < 1.2 x 10“® G [H]. 
Constrained simulations of IGMF [40] favour the val¬ 
ues of R ~ (10“^® — 10“^^) G in voids (see however 
Refs. [¥TI - H5] . advocating somewhat larger values from 
unconstrained simulations and semi-analitical models). 
The angular correlation function of FERMI-LAT pho¬ 
tons points to R ~ 5 X 10“^“^ G jUj. There exist several 
claims of observations of the pair halo around gamma- 
ray sources which suggest IGMF values in the range of 
B ~ 10“^"* G, see e.g. Ref. [35| for a recent one, but these 
analyses are technically involved and require further con¬ 
firmation. 


V. FUTURE TESTS 

While all three methods to distinguish between the two 
scenarios, discussed in Sec. |IV| favour weakly the galac¬ 
tic conversion mechanism, it is clear that more tests are 
required both to confirm the anomalous transparency of 
the Universe and to single out its explanation. To ap¬ 
proach the tests on more solid grounds, future observa¬ 
tions are necessary. Of particular importance are spectral 
and anisotropy studies, for which the following directions 
are especially important: 

• to enlarge the overall statistics of TeV blazars, 
which is best achieved with the coming Gerenkov 
Telescope Array (CTA) |i5] : 

• to study absorption effects in the spectra of 
the most distant blazars, for which one needs 
high-sensitivity observations at energies ~ (10 — 
100) GeV. The sensitivities of both FERMI LAT 
and CTA are insufficient in this energy range; the 


solution may be provided by high-altitude low- 
threshold Cerenkov detectors [47]. Presently, two 
projects of this kind are under consideration, the 
Atmospheric Low Energy Gamma-Ray Observa¬ 
tory (ALEGRO) in Atacama, Chile, and the Elbrus 
Gamma-Ray Observatory (EGO) at the Mount El¬ 
brus, Russia; 

• to move into the strong-absorption energy range 
for bright nearby blazars, which would require 
observations at ~ 100 TeV. The proper instru¬ 
ments for that would be extensive-air-shower de¬ 
tectors, in particular, the Tunka Advanced Instru¬ 
ment for cosmic-ray physics and Gamma-ray As¬ 
tronomy (TAIGA) |35j and the upgraded Carpet 
array at the Baksan Neutrino Observatory |1^ in 
the nearest future. Several years later, the Large 
High Altitude Air Shower Observatory (LHAASO) 
will provide the best sensitivity m- Another rele¬ 
vant future project is HiSCORE El- 

Additional important contributions to the discussion 
are expected from observational constraints on the IGMF 
values. Note that the improvement in the limits on B 
by a factor of two would be sufficient to independently 
exclude the intergalactic scenario. Such an improve¬ 
ment may be achieved after the data of the planned 
all-sky rotation-measure grid [5H [S^, planned for the 
Square Kilometer Array (SKA), will become available. 
Of course, constraints on ga-y-y from laboratory searches 
for the responsible ALP would be crucial, with the most 
sensitive planned instrument being the International Ax- 
ion Observatory (lAXO) [54], for which the entire range 
of photon-ALP couplings suggested in both scenarios we 
discussed is within the discovery range. The upgraded 
Any Light Particle Search (ALPS-II) [5S] experiment will 
approach the interesting range of couplings and probe a 
part of the parameter space for the galactic scenario. 
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